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I.  INTRODUCTION 


He  are  interested  here  in  the  detailed  kinetics  governing  premixed, 
laminar,  one-dimensional  methane/air  flames.  The  earliest  such  scheme'  , 
Introduced  by  Smoot,  Hecker  and  Williams*  and  by  Tsatsaronis,^  Involved  14 
species  and  about  30  reactions.  These  Include  the  reaction  01.0  +  M  ♦  CO  +  H» 
+  H,  which  is  probably  not  an  elementary  reaction  in  atmospheric  pressure 
flames.  Harnatz3  Introduced  a  much  more  complicated  scheme  involving  23 
species  and  SO  reactions.  In  particular,  he  claimed  that  C2  species  were 
Important  for  modeling  rich  or  even  near  stoichiometric  flames.  Dixon-Lewis^ 
attempted  t  >  Identify  the  major  reaction  channels  in  large  mechanisms  such  as 
that  of  Karaatz.  He  compares  two  mechanisms,  one  with  14  species  and  one  with 
18  species  (Including  C «  chemistry).  Both  schemes  were  considered  valid  for 
lean  or  stoichiometric  flames.  Dlxon-Lewis  agrees  with  Warnatz  that  it  is 
necessary  to  Include  C2  chemistry  to  adequately  model  rich  flames. 

In  this  paper  we  introduce  two  new  models.  The  simpler  mechanism 
(14  species,  39  reactions)  does  not  Include  C2  chemistry,  but  unexpectedly 
still  does  a  reasonable  job  of  reproducing  experimental  data  for  lean  and  rich 
flames.  The  more  complicated  model  (20  species,  63  reactions)  is  also  known 
to  be  only  an  approximation,  since  a  number  of  species  and  reactions  have  been 
omitted.  Nevertheless,  it  reproduces  the  experimental  data  quite  well. 

Comparisons  are  made  among  the  seven  kinetic  models.  We  conclude  that 
there  is  insufficient  data  to  validate  any  particular  mechanism  at  this  time. 


II.  NUMERICAL  PROCEDURE 

The  governing  equations  for  the  chemical  species  and  the  temperature  are 
given  in  Reference  5.  These  equations  are  Integrated  in  time,  using  a  finite 
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Flames,"  18th  International  Combustion  Symposium,  The  Combustion  Institute, 
Pittsburgh,  PA,  pp.  369-384,  1981. 

4G.  Dixon-Lewie,  "Aspects  of  the  Kinetic  Modeling  of  Methane  Oxidation  in 
Flames,"  let  Specialists  Meeting  ( International )  of  the  Combustion  Institute, 
France,  pp.  284-289. 

6 T.P .  Coffee  and  J.M.  Hsimerl,  "Transport  Algorithms  for  Premixed,  Laminar 
Steady-State  Flames."  Combustion  and  Flame.  Vol.  43,  pp.  273-289,  1981. 


element  method,  until  the  steady-state  solution  is  obtained, 
calculations  have  been  done  on  a  Cyber  76. 


All 
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Besides  the  kinetics  information,  thermodynamic  and  transport  data  are 
required.  For  most  species,  the  polynomial  fits  of  Gordon  and  McBride*  are 
used.  For  species  not  included  in  the  above,  polynomial  fits  are  made  to  the 
data  of  Benson.^ 

The  molecular  parameters  a  ,  e  /k,  and  p  are  from  Warnatz.3  The 
polarizabilities  are  from  Reid  and  Sherwood. H  These  are  used  to  compute  the 
individual  species  thermal  conductivities  and  the  binary  diffusion 
coefficients.  A  test  problem  was  also  computed  using  the  transport  parameters 
of  Svehla. 12  The  results  were  virtually  identical  with  the  previous 
calculation  using  the  Warnatz  parameters. 

The  multicomponent  mass  fluxes  (pY^V^)  and  heat  flux  (q)  are  computed 
using  the  methods  of  Reference  5.  All  the  preliminary  work  was  done  using 
Method  V,  the  simplest  procedure.  Convergence  took  from  2  to  12  minutes  of 
computer  time,  depending  on  the  mechanism  and  the  stoichiometry.  In  general, 
very  rich  or  very  lean  flames  (near  the  extinction  point)  took  longer  to 
converge.  Once  a  particular  problem  had  been  solved  using  one  kinetics  set, 
the  steady  state  solution  could  be  used  as  the  start  of  the  time  Integration 
for  the  other  kinetics  sets.  This  led  to  substantial  savings. 


6  T.P.  Coffee  and  J.M.  Heimerl ,  "A  Method  for  Computing  the  Flame  Speed  for  a 
Laminar ,  Premixed,  One  Dimensional  Flame,"  ARBRL-TR-02212,  January  1980 
(AD  A08280S). 

^  J.M.  Heimerl  and  T.P.  Coffee,  " The  Detailed  Modeling  of  Premixed,  laminar 
Steady-State  Hamee.  I.  Ozone, "  Combustion  and  Home.  Vol .  39,  pp.  301-315, 
1980. 

8  T.P.  Coffee,  "A  Computer  Code  for  the  Solution  of  the  Equations  Governing  a 
Laminar,  Premixed,  One-Dimensional  Flam,"  AFBRL -MR-031 65,  April  1982. 

8  S.  Gordon  and  B.J.  MeBride ,  " Computer  Program  for  Calculation  of  Complex 
Chemical  Equilibrium  Compositions,  Rocket  Performance,  Incident  and 
Reflected  Shocks  and  Chapman-Jouguet  Detonations,"  NASA-SP-273,  1971  (1981 
program  version). 

10S.  Benson,  Thermochemioal  Kinetics.  2nd  edition ,  John  Wiley  and  Sons,  NY, 
1976. 

11R.C.  Reed  and  J.K .  Sherwood,  The  Properties  of  Gases  and  Liquids.  2nd 
edition,  McGraw-Hill,  NY,  1966. 
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R.A.  Svehla,  "Estimated  Viscosities  and  Thermal  Conductivities  of  Oases  at 
High  Temperatures,"  NASA  Technical  Report  R-l 32,  Lewis  Research  Center , 
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The  flame  speeds  and  profiles  reported  In  this  paper  were  computed  ^sing 
transport  Method  VI.  This  algorithm  Is  very  accurate,  but  is  computationally 
much  more  complex.  Even  using  the  Method  V  solutions  as  starting  points,  the 
14  species  model  required  an  additional  30  to  60  minutes  of  computer  time. 

The  20  species  model  required  an  additional  2  to  4  hours. 


III.  KINETIC  SCHEME 


The  forward  reactions  used  in  this  paper  are  listed  in  Table  1.  The 
rates  of  the  back  reactions  are  computed  from  the  forward  rate  coefficients 
and  the  equilibrium  constant  for  the  reaction.  Additional  species  and 
reactions  were  considered,  but  they  turned  out  to  have  no  noticeable  effect  on 
the  problems  considered. 


The  first  23  reactions  are  from  Dixon-Lewis.  These  have  been  checked 
extensively  in  earlier  studies  of  H2/O2/N2  flames^  and  CO/H2/O2/N2  fl  mi  3.^ 


Reactions  24  through  31  are  derived  from  a  series  of  shock  tube 
experiments  by  Dean,  et  al.15,16  The  formaldehyde  reaction  rates  (28  though 
31)  are  substantially  lower  than  the  rates  previously  used  in  flame  modeling. 


The  CHO  reactions  (32  through  35)  are  not  well  known.  We  have  more  or 
less  arbitrarily  chosen  the  values  of  Dixon-Lewis. ^ 


The  next  four  reactions  (36  through  39)  are  a  miscellaneous  set  from  a 
kinetics  data  base  compiled  by  Gelinas.^ 


These  39  reactions  form  a  14  species  model.  Methane  is  oxidized  along 
the  standard  pathway  CH^  ♦  CH-j  ♦  CHjO  ♦  CHO  +  CO  ♦  CO2 . 

,14 


For  reaction  40,  we  begin  with  the  high  pressure  rate  k  *  2.8  x  10 
T~0.4<  This  rate  is  used  by  both  Warnatz^  and  Dixon-Lewis.^  It  must  be 
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TABLE  1. 

REACTIONS  IN 

THE  CH4/AIR 

SYSTEM 

Reaction 

A* 

B 

C 

Reference 

1. 

OH  +  H2  -  H20  +  H 

1.17E9** 

1.3 

1825 

4 

2. 

H  +  02  -  OH  +  0 

1.42E14 

0.0 

8250 

4 

3. 

0  +  H2  ~  OH  +  H 

1.80E10 

1.0 

4480 

4 

4. 

H  +  02  +  M’  -  H02  +  M' 

1.03E18 

-0.72 

0 

4 

5. 

H  +  H02  -  OH  +  OH 

1.40E14 

0.0 

540 

4 

6. 

H  +  H02  -  0  +  H20 

1.00E13 

0.0 

540 

4 

7. 

H  +  H07  -  H2  +  02 

1.25E13 

0.0 

0 

4 

8. 

OH  +  H02  -  H20  +  02 

7.50E12 

0.0 

0 

4 

9. 

0  +  H02  ~  OH  +  02 

1.40E13 

0.0 

540 

4 

10. 

0  +  H02  -  OH  +  02 

1.25E12 

0.0 

0 

4 

11. 

H  +  H  +  H2  -  H2  +  H2 

9.20E16 

-0.6 

0 

4 

12. 

H  +  H  +  N2  -  H2  +  N2 

1.00E18 

-1.0 

0 

4 

13. 

H  +  H  +  02  -  H2  +  N2 

1.00E18 

-1.0 

0 

4 

14. 

H  +  H  +  H20  ~  H2  +  02 

6.00E19 

-1.25 

0 

4 

15. 

h  +  h  +  co~h2  +  co 

1.00E18 

-1.0 

0 

4 

16. 

H  +  H  +  C02  -  H2  +  C02 

5.49E20 

-2.0 

0 

4 

17. 

H  +  H  +  CH4  ~  H2  +  CH* 

5.49E20 

-2.0 

0 

4 

18. 

H  +  OH  +  M"  -  H20  +  M" 

1.60E22 

-2.0 

0 

4 

19. 

H  +  0  +  M”  -  OH  +  M" 

6.20E16 

-0.6 

0 

4 

20. 

OH  +  OH  -  0  +  H20 

5.75E12 

0.0 

390 

4 

21. 

OH  +  CO  -  C02  +  H 

1.50E7 

1.3 

-385 

4 

22. 

0  +  CO  +  M*  -  C02  +  M’ 

5.40E15 

0.0 

2300 

4 

*  A  is  in  units  of  cm^/mole-sec  or  cm^/mole^-sec,  k«ATBexp  (-C/T). 
**  1.17E9  -  1.17  x  109. 
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23. 

H  +  CO  +  M'  -  CHO  +  M' 

5.00E14 

0.0 

755 

4 

24. 

CH4  +  0  -  CH3  +  OH 

4.07E14 

0.0 

7040 

15 

25. 

CH4  +  H  -  CH3  +  H2 

7.24E14 

0.0 

7590 

15 

26. 

CH4  +  OH  -  CH3  +  H20 

1.55E6 

2.13 

1230 

15 

27. 

ch4  +  m~ch3+h  +  m 

4.68E17 

0.0 

46910 

15 

28. 

CH3  +  0  -  CH20  +  H 

6.02E13 

0.0 

0 

15 

29. 

CH20  +  0  -  CHO  +  OH 

1.82E13 

0.0 

1550 

15 

30. 

CH20  +  H  -  CHO  +  H2 

3.31E14 

0.0 

5290 

15 

31. 

CH20  +  OH  -  CHO  +  H20 

7.58E12 

0.0 

72 

15 

32. 

CHO  +  02  -  CO  +  H02 

3.00E12 

0.0 

0 

4 

33. 

CHO  +  H  -  CO  +  H2 

4.00E13 

0.0 

0 

4 

34. 

CHO  +  Oil  -  CO  +  H20 

5.00E12 

0.0 

0 

4 

35. 

CHO  +  0  ~  CO  +  OH 

1.00E13 

0.0 

0 

4 

36. 

CH20  +  CH3  -  CHO  +  CH4 

2.23E13 

0.0 

2590 

16 

37. 

CH3  +  OH  -  CH20  +  H2 

3.98E12 

0.0 

0 

16 

38. 

ch3  +  ho2  —  ch4  +  o2 

1.02E12 

0.0 

200 

16 

39. 

CO  +  H02  ~  C02  +  OH 

1.50E14 

0.0 

11900 

16 

40. 

CH3  +  Ql3  —  C2H6 

4.56E37 

-7.65 

4250 

4.17 

41. 

C2H6  +  0  ~  C2H5  +  OH 

2.51E13 

0.0 

3200 

16 

42. 

C2H6  +  H  -  C2H5  +  ^ 

5.00E2 

3.5 

2620 

16 

43. 

C2H6  +  OH  -  C2H5  +  H20 

6.63E13 

0.0 

675 

16 

44. 

c2h5  +  h  -  c2h6 

7.23E13 

0.0 

0 

16 

45. 

C2H5  +  H  -  CH3  +  CH3 

3.73E13 

0.0 

0 

16 

46. 

c2h5  -  C2H4  +  H 

2.29E11 

0.0 

19120 

16 

47. 

c2H5  +  o2  —  c2h4  +  ho2 

1.53E12 

0.0 

2446 

16 

48. 

c2h4  +  0  -  ch2  +  CH20 

2.53E13 

0.0 

2516 

16 

49. 

C2H4  +  OH  -  CH20  +  CH3 

5.00E13 

0.0 

3020 

16 

13 


2.53E13 


0.0 


2516 


16 


50. 

C2H4  +  0  -  C2H3  +  OH 

2.53E13 

0.0 

2516 

16 

51. 

c2h4  +  o2  #  c2h3  +  ho2 

1.3  3E 1 5 

0.0 

27680 

16 

52. 

C2H4  +  H  -  C2H3  +  Hj 

2.00E15 

0.0 

10000 

16 

53. 

C2H4  +  oh  -  c2h3  +  h2o 

4.40E14 

0.0 

3270 

16 

54. 

c2h3  +  m  -  c2h2  +  h  +  m 

3.01E16 

0.0 

20380 

16 

55. 

c2H3  +  02  -  C2H2  +  H02 

1.57E13 

0.0 

5030 

16 

56. 

C2H3  +  H  -  C2H2  +  Hj 

7.53E13 

0.0 

0 

16 

57. 

c2h3  +  oh  -  c2h2  +  h2o 

I.00E13 

0.0 

0 

16 

58. 

C2H2  +  OH  -  CH3  +  CO 

5.48E13 

0.0 

6890 

16 

59. 

CH3  +  h  **  ch2  +  h2 

2.00E11 

0.7 

-1500 

16 

60. 

CH3  +  OH  -  CH2  +  H20 

6.00E10 

0.7 

1010 

16 

61. 

CH2  +  02  -  CHO  +  OH 

1.00E14 

0.0 

1860 

16 

62. 

ch2  o2  “  ch2o  +  0 

1.00E14 

0.0 

1860 

16 

63. 

ch2  +  02  •  C02  +  H2 

1.00E14 

0.0 

1860 

16 

(Ml 

»  total  concentration 

l  M  *  J 

-  [H2]  +  0.74  [CO]  +  l.i 

o 

+ 

CM 

8 

.35  [021  +  6.5 

[H201  +  0.44 

[n21 

[M"l 

*  Ih2)  +  [CO]  +  [co2l  + 

[02]  +  5.0 

[h2o!  +  [n2I 
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further  modified  by  the  fall  off  method  of  Luther  and  Troe.  °  The  parameters 
in  Table  1  are  a  least  squares  fit  of  the  modified  values  for  a  pressure  of 
one  atmosphere.  This  is  the  primary  channel  for  the  formation  of  C2 
hydrocarbons.  The  C2  reactions  are  40  through  58.^  The  primary  pathway  is 


C2h6  *  c2h5  *  c2h4  *  c2h3  *  C2H2 


CH, 


CO 


There  is  al-.  o  a  side  channel  through  reactions  48  and  49, 


02^  ^  CH2O. 

The  last  five  reactions  (59  through  63)  consist  of  CH2  chemistry.1?  The 
pathway  is 


ch3  ♦  ch2 


CHO 


CO. 


This  sequence  has  not  been  included  in  previous  CH^/air  flame  models. 


IV.  COMPARISON  WITH  EXPERIMENTAL  DATA 


Here  we  compare  the  seven  models  discussed  above  with  experimental 
data.  All  flames  considered  in  this  section  are  atmospheric  pressure, 
unbounded,  methane/air  flames  with  an  initial  temperature  of  298K. 

Many  measurements  have  been  made  on  the  burning  velocity  of  methane/alr 
flames.  Andrews  and  Bradley19  have  critically  reviewed  the  different 
experimental  techniques.  Following  their  recommendations,  we  have  used  the 


18K.  Luther  and  J.  Troe,  ” Weak  Collision  Effects  in  Dissociation  Reactions  at 
High  Temperatures,"  17th  International  Combustion  Symposim,  The  Combustion 
Institute,  Pittsburgh,  PA,  pp.  535-542 ,  1979. 

19G.E.  Andrews  and  D.  Bradley,  " Determination  of  Burning  Velocities:  A 
Critical  Review,"  CambrntiaiLand  Flame,  Vol .  18,  pp.  133-153,  1972. 


results  of  four  experiments. 20-23  uhiie  there  is  scatter  in  the  data,  the 
results  are  generally  consistent. 


The  node Is  of  Smoot,  Hecker  and  Williams  (SHW14)  and  Dixon-Lewis  <DL14 
and  DL18)  are  accurate  for  lean  to  slightly  rich  flames.  Tsatsaronis  [TS14] 
modified  the  model  of  Smoot,  et  al,  so  that  it  is  accurate  for  rich  flames, 
but  the  calculated  burning  velocities  are  low  for  stoichiometric  to  slightly 
lean  flames.  The  more  complicated  model  of  Warnatz  (W23)  is  reasonably 
accurate  over  the  entire  range,  although  the  flame  speeds  are  low  for  rich 
flames. 

Figure  1  shows  the  results  for  the  two  new  models  introduced  in  this 
paper  (CF14  and  CF20),  the  model  of  Warnatz,  and  the  experimental  burning 
velocities.  The  results  are  given  in  terms  of  the  equivalence  ratio 
$  *  2  XCHa/X02»  where  X  represents  the  initial  mole  fraction  of  the  species. 
All  experimental  values  have  been  corrected  to  298  K.  Both  new  models  are 
accurate  over  the  entire  range,  although  slightly  low  for  very  lean  or  very 
rich  flames.  The  values  calculated  here  using  the  Warnatz  model  (W23)  for 
rich  flames  »re  slightly  lower  than  those  reported  by  Warnatz. ^  We  do  not 
know  the  reason  for  the  discrepancy. 

Laser  methods  have  recently  been  used  to  measure  species  and  temperature 
profiles.  Bechtel  and  Teets^  measured  OH  profiles  using  laser-induced 
fluorescence  for  $  -  0.86,  1.00,  and  1.25.  In  a  later  paper,  Bechtel, 
et  al ,  measured  H2,  CO,  CH^,  O2,  CO2 ,  H2O,  and  temperature  profiles  for  the 
same  three  cases. 

These  profiles  were  compared  with  calculations  using  the  SHW14  model. 
Agreement  was  very  good  except  for  the  H2  profile.  The  measured  H2 
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G.E.  Andrew  and  D.  Bradley,  "Determination  of  Burning  Velocity  by  Double 
Ignition  in  a  Cloeed  Vessel,"  Combustion  and  Flame,  Vol.  20,  pp.  77-89, 

1973. 

pi 

R.  Gunther  and  G.  Janish,  "Measurements  of  Burning  Velocity  in  a  Flat  Flame 
Front,"  Combustion  and  Flame .  Vol.  19,  pp.  49-53,  1972. 
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R.  Lindow,  " Eine  verbeseerte  Brennermethode  sur  Bestimmg  der  lamina^en 
flajmengeschwindigkeiten  von  Brenngae/Luft-Gemiechen , "  Bremetaff  Warme 
Kraft.  Vol.  20,  pp.  8-14,  1968. 
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S. B.  Reed,  J.  Mineur,  and  J.P.  McNaughton,  "The  Effect  on  the  Burning 
Velocity  of  Methane  of  Vitiation  of  Combustion  Air,"  J.  Inst.  Fuel,  Vol.  44, 
pp.  149-155,  1971. 

Bechtel  and  R.E.  Teets,  "Hydroxyl  and  Its  Concentration  Profile  in 
Methane-Air  Flame."  Applied  Ovtics.  Vol ,  18,  pp.  4138-4144,  1979. 

^sd.H.  Bechtel,  R.J.  Flint,  C.J.  Dasch ,  and  D.A.  Weinberger,  "Atmospheric 
Pressure  Premixed  Hydrocarbon-Air  Flams:  Theory  and  Experiment, " 
Combustion  and  Flame.  Vol.  42,  pp.  197-213,  1981. 


concentration  was  such  greater  than  the  calculated  H2  concentration  in  the 
leading  edge  of  the  flame. 


One  explanation  considered  by  the  authors  is  an  inadequate  treatment  of 
the  diffusion  of  Hj  in  their  model.  This  appears  to  be  the  case.  When  the 
SHW14  model  was  run  using  our  code,  agreement  for  the  Hj  profile  was  very 
good . 


Calculations  were  made  for  each  of  the  three  cases  above  for  all  seven 
kinetic  models.  Agreement  was  excellent  for  all  models  for  the  CO2,  02,  CH^, 
H2O,  and  temperature  profiles.  The  primary  difference  was  that  the 
experimental  temperature  dropped  in  the  post-flame  region,  while  the  model 
temperatures  increased.  This  can  be  explained  by  heat  lost  to  the  wall. 

Agreement  for  the  H2  and  CO  profiles  was  very  good,  but  with  some  scatter 
in  the  calculated  peak  heights. 

Agreement  for  the  OH  profile  was  good  for  the  lean  and  stoichiometric 
cases.  For  the  rich  flame,  all  of  the  models  predicted  a  higher  and  earlier 
peak  than  that  given  by  the  experiment.  This  suggests  either  a  problem  with 
the  experiment  or  some  basic  inadequacy  in  the  standard  models  for  rich  flames 
or  both. 

Figures  2-5  show  the  experimental  profiles  for  a  stoichiometric  flame 
compared  to  the  3  models  Including  C2  chemistry.  The  other  models  are  quite 
similar. 


V.  SENSITIVITY  ANALYSIS 


A  sensitivity  analysis  is  useful  in  trying  to  understand  a  complex  system 
like  the  CH^ 'air  kinetics.  A  procedure  for  finding  sensitivity  coefficients 
for  flames  vas  developed  in  Reference  26. 


This  analysis  was  done  for  a  stoichiometric  CH^/air  flame  using  the  CF20 
model.  The  larger  logarithmic  sensitivity  coefficients  Sg*  for  the  flame 
speed  are  given  in  Table  2.  These  are  defined  such  that  If  the  rate  of  the 
reaction  is  changed  by  a  small  factor  a,  the  flame  speed  will  change  by  a  to 


the  power,  i.e.,  v  =  vn  a 

0.5.2<F  0 


Analytically,  the  sian  of  Sg1  must  equal 


The  flane  is  most  sensitive  to  the  chain  branching  reaction  2  IH  + 

(>2  +  OH  +  Hj.  Radicals  must  be  produced  for  the  flame  to  propagate,  and 
reaction  2  is  the  principal  bottleneck  for  reaction  production. 

The  fl;  me  is  much  less  sensitive  to  the  corresponding  chain  branching 
reaction  3  1 l  +  Hj  OH  +  H] .  Since  the  H2  concentration  is  small,  this 


Ooffei  and  J.M.  Heimerl ,  "Sensitivity  Analysis  for  Premixed  Laminar , 
Steady-State  Hams,"  Combustion  and  Plante ,  Vol,  SO,  pp.  323-340,  1983 . 


Figure  2.  Stoichiometric  CH^/Air  Flame.  Model  Values  -  CF20  (Line),  W23  (Dot),  and  DL18  (Dash) 

Bechtel  Experimental  Points.  T  and  OH  Profiles. 


Figure  3.  Stoichiometric  CH^/Air  Flame.  Model  Values  -  CF20  (Line),  W23  (Dot),  and  DL18  (Dash) 

Bechtel  Experimental  Points.  H2  and  CO  Profiles. 


Figure  5.  Stoichiometric  CH^/Air  Flame.  Model  Va 

Bechtel  Experimental  Points 


TABLE  2.  LOGARITHMIC  SENSITIVITY  COEFFICIENTS  FOR  A 
STOICHIOMETRIC  FLAME  (CF20  KINETICS  SCHEME) 


Reaction 

SEJ 

Reaction 

sEJ 

1 

.06 

40 

1 

e 

O 

2 

.45 

59 

.09 

3 

.04 

60 

.03 

4 

-.15 

61 

.03 

18 

-.05 

62 

.03 

21 

.16 

63 

-.06 

23 

.07 

25 

-.06 

27 

-.06 

28 

.03 

33 

1 

• 

O 

reaction  is  relatively  unimportant.  Most  of  the  H  atoms  are  generated  by 
other  reactions. 

Much  of  the  sensitivity  analysis  can  be  explained  in  terms  of  a  shortage 
of  H  atoms.  Reactions  that  consume  H  atoms  [4,  18,  25,  27,  33]  tend  to  slow 
down  the  flame.  Reactions  that  produce  H  atoms  [1,  3,  21,  23,  28]  increase 
the  flame  speed.  The  flame  is  much  more  sensitive  to  the  radicals  than  to  the 
rate  at  which  the  intermediates  [CH3,  C^O,  CHO,  00}  are  formed. 

Reaction  40  [CH-j  +  CH^  ♦  02^]  slows  down  the  flame  because  it  leads  to 
the  slower  C~  pathway.  However,  only  about  10X  of  the  CH3  formed  goes  to 
C2H*.  All  trie  other  C2  chemistry  rates  have  very  small  sensitivity 
coefficients.  For  other  kinetic  schemes,  more  CHj  may  be  converted  to  C2H^, 
and  the  C2  reactions  would  be  more  important. 

Finally,  reactions  59  fCH3  +  H  ♦  CH2  +  H2]  and  60  (CH3  +  OH  ♦  CH2  +  H20] 
have  positive  sensitivity  coefficients,  even  though  radicals  are  consumed. 

This  is  because  CH2  reacts  rapidly  with  02  to  produce  radicals  through 
reactions  61  and  62*  The  corresponding  reaction  63,  which  does  not  produce 
radicals,  has  a  negative  sensitivity  coefficient. 


In  summary,  the  flame  speed  Is  moat  sensitive  to  the  production  of 
radicals.  The  kinetics  Is  the  most  Important  for  determining  the 

flame  speed.  The  next  most  Important  group  is  the  00  to  CO2  chemistry 
(reactions  21  and  22).  In  contrast,  the  flame  Is  fairly  Insensitive  to  the 
specific  hydrocarbon  reactions. 

The  Cj  chemistry  (reactions  40  to  58)  tends  to  slow  down  the  flame,  while 
the  CHj  to  CH2  chemistry  (reactions  59  to  63)  increases  the  burning 
velocity.  The  latter  reactions  are  not  standard  for  CH^  chemistry  models. 

But  if  these  reactions  are  omitted,  the  burning  velocity  is  too  low, 
especially  for  rich  flames.  So  these  reactions  are  important  for  the  CF20 
mechanism.  Whether  these  reactions  are  actually  important  in  CH^  combustion 
is  yet  to  be  determined. 

Sensitivity  coefficients  for  the  species  and  temperature  profiles  were 
also  computed  (as  a  function  of  position).  The  major  species  (CH4,  O2,  C02» 
H2O)  were  the  least  sensitive  to  changes  in  the  rate  coefficients.  Next  comes 
the  intermediate  CO  and  H2  and  the  OH  radical.  The  other  radicals  (H,  0,  HO2) 
were  somewhat  more  sensitive.  The  intermediates  along  the  main  oxidation 
pathway  (CH-j,  CH2O,  CHO)  were  quite  sensitive.  The  species  most  affected  by 
rate  changes  were  the  C2  species  and  CH2« 

Moreover,  the  species  CH4,  O2,  CO2,  H2O,  CO,  H2,  OH,  H,  and  0  are  most 
sensitive  to  the  rate  constants  2  and  4.  It  is  the  species  along  the  main 
oxidation  pathway  (CH-j,  CT^O,  and  CHO)  that  are  sensitive  to  the  hydrocarbon 
reactions. 

This  explains  why  all  the  models  are  in  very  good  agreement  with  the 
profiles  measured  by  Bechtel.  The  profiles  measured  are  the  least  sensitive 
profiles,  and  in  any  case,  are  primarily  sensitive  to  changes  in  the  H2/O2 
subset.  A  wide  range  of  hydrocarbon  reactions  and  rates  can  be  assumed,  as 
long  as  an  appropriate  pathway  exists  from  CH^  to  CO  to  C02  and  the  flame 
speed  is  reasonable. 


VI.  DISCUSSION 

There  is  substantial  agreement  among  the  models  for  the  profiles  measured 
by  Bechtel,  et  al.  It  is  also  useful  to  see  for  what  profiles  the  models 
differ.  Experimental  measurements  of  such  profiles  could  help  determine 
which,  if  any,  of  the  models  are  correct. 

To  help  analyze  the  results,  a  screening  analysis  is  performed. ^  That 
is,  the  rate  of  production  and  loss  for  each  species  as  a  function  of  position 
is  partitioned  according  to  the  contribution  of  each  reaction.  This  shows  the 
pathways  by  which  the  various  species  are  produced  and  consumed  in  the  given 
network.  Screening  analysis  complements  a  sensitivity  analysis,  which  shows 
how  rate  changes  can  effect  the  given  network. 

First,  we  consider  the  importance  of  the  C£  and  CH2  chemistry.  Table  3 
gives  the  ratios  of  the  peak  mole  fractions  of  the  CF20  and  CF14  models  as  a 
function  of  stoichiometry,  as  well  as  the  flame  speeds  S.  The  major  species 
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TABLE  3. 

RATIO  OF 

CF20  PEAK 

HEIGHTS  TO  CF14  PEAK  HEIGHTS 
EQUIVALENCE  RATIO 

AS  A  FUNCTION  OF 

0.65 

0.86 

1.00 

1.25 

1.40 

1.50 

H 

.86 

.94 

.95 

.91 

1.15 

1.33 

OH 

.98 

1.00 

1.00 

.96 

1.14 

1.06 

0 

.92 

.96 

.96 

.92 

1.86 

1.87 

CM 

s 

1.57 

1.41 

1.42 

1.70 

2.20 

1.81 

ch3 

.53 

.52 

.49 

.42 

.41 

.31 

ch2o 

.89 

.99 

1.07 

1.23 

1.51 

1.59 

CHO 

.80 

.93 

1.00 

1.09 

1.39 

1.09 

H2 

.94 

1.02 

1.04 

1.05 

1.02 

1.04 

00 

.88 

.92 

.92 

.91 

.94 

1.01 

s 

.95 

1.07 

1.12 

1.17 

1.26 

1.07 

(CH^,  O2,  C02,  H2O)  and  the  temperature  are  not  included  as  the  profiles  are 
virtually  identical.  The  exception  is  the  rich  flames  ($  *  1.40  or  $  «  1.50). 
where  the  CO2  profiles  show  slight  differences. 

For  lean  to  slightly  rich  flames,  only  the  CH3  and  H02  profiles  show 
large  differences.  The  CF20  model  has  less  CH3  since  the  CH^  c2H6  and 
CHj  +  CH2  reactions  deplete  this  species.  The  C2  chemistry  then  leads  to  more 
HO2  through  the  reactions  48  {C2H5  +  O2  ♦  C2H4  +  H021  and  56  [C2H3  +  02  *  C2H2 

+  ho2] . 

The  similarities  between  the  two  models  are  due  to  competing  effects. 

The  C2  chemistry  tends  to  slow  down  the  flame.  The  bulk  of  the  C2H^  formed 
slowly  oxidizes  to  C-^,  and  then  goes  to  00  through  reaction  58  [C2H2  + 

OH  ♦  CH-j  +  CO] .  The  concentrations  of  the  radicals  and  the  Intermediates  CH20 
and  CHO  are  all  lowered. 

The  CH2  chemistry  tends  to  increase  the  flame  speed.  Once  the  CH2  forms, 
it  oxidizes  fairly  rapidly  by  the  reactions  61  [ CH2  +  02  ♦  CHO  +  OH],  62  [CHo 
+  02  ♦  CH20  +  0]  and  63  (CH2  +  02  ♦  COo  +  H2I.  The  intermediates  CH20  and  CHO 
are  Increased.  The  radical  concentrations  are  not  greatly  affected,  although 
the  0  atom  peak  does  increase. 

For  the  rich  flames,  more  of  the  CHj  goes  to  C2H6  (up  to  20%)  and  also  to 
CH2  (up  to  35%).  This  leads  to  slightly  higher  flame  speeds  for  the  CF20 
model.  The  CH3  peak  becomes  much  smaller.  The  0  and  CH20  concentrations  are 
substantially  higher,  due  to  reaction  62  [CH2  +  02  ♦  +  0]. 
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The  Hj  profile  shows  relatively  minor  differences.  For  the  CF14  model, 
the  H2  is  primarily  produced  through  reactions  25  [ CH^  +  H  ♦  CH-j  +  H2 ]  and  30 
lffl20  +  H  ♦  CHO  +  H2).  For  the  CF20  model,  less  H2  is  produced  through  these 
reactions,  but  this  is  compensated  for  by  reactions  59  [CH^  +  H  ♦  CH2  +  H2] 
and  42  (C2H^  +  H  ♦  C2Hg  +  H2J .  The  overall  rate  of  production  through  the 
various  pathways  is  about  tne  same.  As  long  as  the  flame  speeds  are  similar, 
this  will  be  the  case. 

Similarly,  a  decrease  in  the  CO  production  reaction  23  [CHO  +  M'  +  H  +  03 
+  M'  ]  is  compensated  for  by  reaction  58  [  C2H2  +  OH  -*■  CHj  +  00] . 

Next,  we  consider  comparisons  among  the  three  models  with  chemistry. 

To  conserve  computer  time,  the  solutions  are  generated  using  the  simplest 
transport  algorithm.  Method  v.^  Tables  4  and  5  show  the  ratios  of  the  peak 
mole  fractions  and  flame  speeds.  Again,  the  major  species  and  temperati  re 
profiles  are  virtually  identical. 

The  flame  speeds  generated  by  the  DL18  model  agree  with  those  reported  by 
Dixon-Lewis.^  The  flame  speeds  generated  by  the  W23  model  are  lower  than 
those  reported  by  Warnatz  (5%-15%). 

For  the  DL18  model,  results  are  not  given  for  $  =1.50.  No  steady  state 
solution  could  be  generated,  as  the  computed  flame  always  went  out. 

For  lean  to  slightly  rich  flames,  all  three  models  agree  for  the 
temperature  and  major  species  profiles.  Lower  flame  speeds  are  predicted 
for  <j>  *  1.40  (DL18  and  W23  models)  and  for  <)>  =  1.50  (W23  model).  At  these 
lower  speeds,  the  flame  front  is  spread  out  and  noticeable  differences  appear 
in  the  profiLes. 

The  H,  OH,  and  0  profiles  are  similar  for  lean  to  stoichiometric 
flames.  For  rich  flames,  where  the  models  predict  lower  flame  speeds,  the 
radical  concentrations  tend  to  be  lower  (up  to  a  factor  of  2) .  This  is 
reasonable  since  large  amounts  of  radicals  lead  to  rapid  combustion. 

There  are  large  differences  in  the  intermediates  CH^,  CH20,  and  CHo.  CHO 
is  especially  variable  because  of  its  low  concentration.  The  C2  species  show 
relatively  low  concentrations  and  large  variation.  But  the  H2  and  00 
profiles,  which  have  larger  concentrations,  remain  similar  throughout. 


TABLE  4. 

RATIO  OF  DL18 

PEAK  HEIGHTS  TO  CF20  PEAK 
EQUIVALENCE  RATIO 

HEIGHTS  AS  A 

FUNCTION  OF 

0.65 

0.86 

1.00 

1.25 

1.40 

H 

.93 

.98 

.99 

1.01 

38 

OH 

.98 

1.00 

1.00 

1.01 

.50 

0 

.96 

1.00 

1.00 

.91 

.23 

ho2 

.54 

.70 

.79 

.78 

.40 

C2Hf) 

3.35 

3.53 

3.31 

2.26 

1  .  10 

C2H3 

.23 

.51 

.77 

1.08 

1.21 

c2h4 

.84 

1.71 

2.80 

3.62 

2.30 

ch3 

1.35 

1.53 

1.62 

1.74 

1.49 

ch2o 

.62 

.64 

.69 

.97 

1.35 

CHO 

2.16 

2.24 

2.11 

1.95 

.81 

«2. 

.79 

.82 

.83 

.87 

.88 

CO 

1.15 

1.14 

1.12 

1.10 

.97 

s 

1.00 

1.00 

.96 

.88 

.50 

VII.  LOW  PRESSURE  FLAMES 


There  have  been  several  experimental  studies  of  low  pressure  flam.  using 
mass  spectroscopy.  We  consider  one  such  experiment  due  to  Fristrom,  e-  al,  on 
a  7.853:  (H4/O2  flame  at  0.05  atm  pressure.  The  temperature,  CH^,  0~,  >2 . 

H20.  H2,  00,  and  CH2O  profiles  were  measured.  The  flame  was  stabilize  on  a 
circular  screen  burner.  The  experimental  burning  velocity  (referenced  o 
298K,  0.05  atm)  was  69  cm/sec. 

The  flame  was  modeled  as  an  unbounded  flame  with  the  temperature  of  the 
unburned  mixture  as  298K.  Computations  were  made  using  three  kinetics  schemes 
(CF20,  W23,  and  DL18)  and  the  more  accurate  transport  Method  VI.  Sine  •  the 
flame  is  very  lean,  the  C2  chemistry  has  only  minor  effects. 

The  temperature  profiles  were  in  very  good  agreement  through  the  flame 
front.  In  the  post  flame  region  the  experimental  temperature  measurements 
were  high.  For  the  major  species  CH^,  O2,  OO2 ,  and  H2O,  all  three  models  were 
in  very  good  agreement  with  the  experimental  results.  There  were  noticeable 
differences  for  the  other  species  peak  heights  and  the  burning  velocities  (see 
Table  6).  Dixon-Lewis^  tried  adjusting  the  rate  constants  to  match  tl 
profiles.  He  was  unable  to  do  so  and  still  match  the  atmospheric  pressure 
experimental  results. 

Fristrom,  et  al,  noted  that  recombination  in  the  sampling  probe  could 
distort  the  trace  species.  To  test  this  possibility,  an  ODE  code  was  used. 

The  initial  temperature  and  concentrations  were  taken  as  the  values  at  the 
peak  OUO  concentration,  as  determined  by  the  CF20  model.  A  time  integration 
was  performed  until  the  CH^O  concentration  matched  the  experimental  value, 
holding  the  temperature  ana  pressure  constant.  The  results  are  given  in  Table 
4.  The  values  for  Hj  and  00  also  closely  match  the  experimental  values.  The 
major  species  concentrations  also  change,  but  since  their  concentrations  are 
larger  to  begin  with,  this  does  not  affect  their  agreement  with  the 
experimental  profiles. 


TABLE  6.  A  COMPARISON  OF  BURNING  VELOCITIES  AND  PEAK  MASS  FRACTIONS 


FOR  THE  LOW 

PRESSURE  FLAME 

Experiment 

CF20 

W23 

DL18 

Burning  Velocity 

69 

55.7 

62.1 

39.5 

CH,0 

1.5  E-3 

3.1  E-3 

1.1  E-3 

1.2  E-3 

2.2  E-4 

2.3  E-4 

1.3  E-4 

9.0  E-5 

c6 

4.1  B-2 

3.4  B-2 

3.6  E-2 

3.2  E-2 
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TABLE  7.  RESULTS  OF  TOE  TIME  INTEGRATION  OF  TOE  CF20  MODEL  PROFILES. 
INITIAL  CONDITIONS  ARE  TAKEN  AT  TOE  CH20  PEAK 


ch2o 

& 


t  ■  o 
3.1  E-3 
2.3  E-4 
3.0  E-2 


t 


g3.2  E-4 
1.5  E-3 
2.2  E-4 
3.8  E-2 


sec 


So  it  is  possible  that  the  differences  between  the  CF20  model  and  the 
experiment  are  due  to  recombination  within  the  probe.  On  the  other  hand, 
there  is  no  obvious  way  recombination  could  lead  to  agreement  between  either 
the  W23  or  the  DL18  scheme  and  the  experiment.  However,  reactions  on  the  wall 
of  the  probe  and  the  cooling  due  to  the  probe  have  been  ignored  in  this  simple 
model.  These  assumptions  cannot  be  justified,  and  the  above  explanation  must 
be  considered  tentative. 


VIII.  CONCLUSIONS 

Results  for  seven  kinetic  schemes  have  been  compared  with  a  range  of  data 
for  atmospheric  pressure  methane/air  flames.  Three  of  the  models  (the  two 
introduced  here  and  the  one  due  to  Warnatz)  are  reasonably  accurate  over  the 
entire  range. 

None  of  the  models  match  the  data  on  low  pressure  flame  obtained  by  mass 
spectroscopy.  However,  if  chemistry  is  assumed  to  continue  for  a  short  time 
in  the  probe,  the  predictions  of  the  two  new  models  are  consistent  with  the 
experimental  data. 

None  of  the  models  can  be  considered  to  be  completely  validated.  In 
particular,  the  CF14  models  ignores  C2  chemistry,  although  it  is  known  that  C2 
species  are  formed.  However,  methane  is  oxidized  in  a  series  of  steps.  As 
long  as  the  rate  of  the  overall  pathway  is  accurate,  the  flame  speed,  the 
temperature  profiles,  and  the  major  species  profiles  will  be  accurate.  But 
this  does  not  guarantee  that  the  intermediates  CH3,  CH20,  CHO,  CH2,  and  the  C2 
species  will  be  correct. 

A  measurement  of  the  CH3  profile  would  be  especially  useful.  This  would 
help  determine  the  branching  ratio  between  CH3  -*■  CH20  and  CH3  ♦  C2H6. 

So  there  are  several  models  that  reproduce  well  the  main  features  of 
atmospheric  flame.  Most  of  the  models  are  accurate  for  lean  to  slightly  rich 
flames.  But  the  details  of  the  intermediate  reactions  in  the  oxidation  of  CH^ 
are  still  not  well  known. 
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